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ABSTRACT 

S" 

°\ 

We present an efficient method for computing lightcurves of an elliptical source 

which is microlensed by a point mass. The amplification of an extended source involves 

a two-dimensional integral over its surface brightness distribution. We show that for 

■^ a general surface brightness profile with an elliptical symmetry, this integral can be 

reduced to one dimension. We derive analytical results for the entire lightcurve in 

the limit of low (e.g. planetary) mass lenses, and for the wings of all microlensing 

lightcurves in general. In both cases the lightcurve carries information about deviations 

of the source from elliptical symmetry, e.g. due to spots. The method is used to find 

the amplification of a circular red giant photosphere and an inclined accretion disk. 

We demonstrate that microlensing of an emission line from a disk can be used to infer 

the disk velocity structure and surface brightness profile. 
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1. Introduction 

The sample of all known gravitational lenses is by now dominated by stellar microlenses in 
our galaxy (Alcock et al. 1995, 1996; Udalski et al. 1994; Ansari et al. 1996; Alard 1996). Even 
galactic macro-lenses of quasars are often accompanied by microlensing events due to stars (Irwin 
et al. 1989; Schild and Smith 1991; Schneider, Ehlers, & Falco 1992 and references therein), 
and thus emphasize that microlensing is a common phenomenon. Based on the fact that a solar 
mass lens at a distance al offers an angular resolution ~ (<i/10Gpc)^ 1 ' 2 micro-arcsec, microlensing 
has long been recognized as a potential tool for resolving the surface brightness distribution of 
compact sources. 

Nevertheless, most stellar sources are smaller than the effective size of their lenses. The finite 
extent of a stellar photosphere (~ Rq ) is still much smaller than the Einstein radius of a solar mass 
lens [~ 10 3 i? Q x (Mi ens /M Q ) 1 / 2 (o!/10kpc) 1 / 2 , where Mi ens is the lens mass], and so the lightcurve 
of galactic microlenses typically carries little information about the brightness distribution of the 
source (Paczyhski 1986; Udalski et al. 1994). The point source approximation breaks down only 



when the projected impact parameter is comparable to the finite size of the source. Although such 
circumstances are realized in a minority of all microlensing events, they carry a large scientific 
payoff. A spectroscopic identification of the source size and distance in such a case, could break 
the parameter degeneracy of the microlensing lightcurve and reveal the Einstein angle of the lens 
(Loeb & Sasselov 1995; Gould h Welch 1996). The appearance of extended source signatures 
should be common when giant stars are microlensed by planets (Loeb & Sasselov 1995; Bennett 
& Rhie 1996). By now, there are some examples of lightcurves with unusually high amplification 
(i.e. small impact parameter) or short duration (i.e. low lens masses) in which there are hints of 
finite source signatures (Mao et al. 1994; Lennon et al. 1996; Bennett et al. 1996; Pratt 1996; 
Alcock et al. 1997). Searches for microlensing by planets (Albrow et al. 1996; Pratt et al. 1996) 
are likely to add more examples to this class of events in the future. Spectroscopic monitoring of 
finite source events could ultimately test the theory of stellar photospheres, especially for giant 
stars (Loeb & Sasselov 1995; Sasselov 1996). 

An analogous ordering of scales appears in stellar microlensing of the continuum emission 
from quasars. For a star of mass Mi cns , the Einstein radius at a cosmological distance is of order 
~ 5 x 10 16 (Miens/A^,) 1 ' 2 cm. In comparison, the optical continuum emission in quasars is believed 
to originate from a compact accretion disk. The UV bump observed in quasar spectra is often 
interpreted as thermal emission at a disk temperature T<ji s k ~ 10 5 Ts K (e.g. Laor 1990), and so 
the scale of the disk emission region must be ~ 10 15 T 5 ~ L 4 g cm, where L46 is the corresponding 
luminosity of the quasar in units of 10 46 erg s _1 . Thus, for lens masses Mi ens S> 1O~ 3 M0, the 
continuum source is much smaller than the projected Einstein radius of the lens. This expectation 
is indeed confirmed in the lenses of Q2237+0305 (Irwin et al. 1989) and 0957+561 (Schild & 
Smith 1991), where variability due to microlensing has been detected (Wambsganss, Paczyhski, &: 
Schneider 1990; Rauch & Blandford 1991; Racine 1991; Gould & Miralda-Escude 1996). 

The above ordering of scales is reversed for the broad emission line regions of quasars. 
Reverberation studies of the time lag between variations in the continuum and the line emission 
in active galactic nuclei (AGN) indicate that the broad emission lines originate at a distance 
of ~ 3 x 10 17 L 4 g cm (Peterson 1993; Maoz 1996). This implies that the broad lines could be 
amplified through stellar microlensing only for low-luminosity AGN (L^q <J 10~ 2 ). Microlensing 
of Seyfert galaxies offers a unique opportunity for resolving the unknown structure and velocity 
distribution of the broad line region in AGN. 

For all of the above applications, it is important to find the amplification in a variety of 
source and lens geometries. Evaluating the net amplification of an extended source involves a 
two-dimensional integral over its surface brightness distribution. A direct calculation of this 
integral raises some numerical problems (due to the divergence of the amplification at the lens 
position) and is computationally time consuming. Witt & Mao (1994) expressed the amplification 
of a circular source with uniform brightness in terms of elliptic integrals. Alternatively, the 
amplification of uniform sources can be computed in the image plane, avoiding the divergence. 
Realistic sources, however, have a non-uniform surface brightness distribution (e.g. due to 



limb-darkening in stellar atmospheres) that also depends on the observed wavelength (Loeb 
& Sasselov 1995; Gould & Welch 1996; Sasselov 1996). Witt & Mao (1994) proposed to sum 
over rings using the derivative of the brightness distribution as the appropriate weight for their 
contributions to the overall amplification. As a result, they traded one dimension in their surface 
integration for the evaluation of their elliptic integral solution. Other methods, such as image 
plane integration, require approximating nonuniform sources as a superposition of uniform sources. 
In this work, we propose another, more efficient approach to microlensing of nonuniform sources. 
We show how to analytically reduce the surface integral to one dimension for elliptically symmetric 
brightness profiles. This generalization can therefore be used to calculate microlensing of inclined 
disks, in addition to its straightforward application to stellar photospheres (Heyrovsky, Loeb, & 
Sasselov 1996). 

Microlensing is also a valuable tool for studying sources with organized internal motions, 
such as rotating disks. If the gas in the disk emits a spectral line, then the variation of the line 
profile during a microlensing event would contain important information about the trajectory of 
the lens and the inclination angle, velocity structure, and brightness profile of the disk. As the 
lens moves across the disk, it amplifies the intensity at the Doppler shifted wavelength associated 
with the rotation velocity behind the lens. By examining the time dependence of the line profile 
it is possible, in principle, to deconvolve the radial brightness profile of the disk down to the 
innermost radius sampled by the lens trajectory. Examples for relevant sources include the iron 
Ka fluorescence line in AGN (Tanaka et al. 1995; Fabian et al. 1995; Iwasawa et al. 1996), 
the maser lines emitted from compact disks in galactic nuclei (e.g., Miyoshi et al. 1995), and 
the line emission from stellar accretion disks (Warner 1995, and references therein). A related 
effect has previously been discussed in the context of the broad emission lines in AGN (Nemiroff 
1988; Wambsganss 1990; Schneider & Wambsganss 1990), although these lines are not believed to 
originate from a cold Keplerian disk because they do not posses a double-peaked profile. 

In this work, we restrict our attention to microlensing by an isolated point mass lens, and 
avoid more complex lensing environments (see, e.g. Jaroszyhski, Wambsganss, & Paczyhski 1992; 
Wambsganss & Kundic 1995). 

The outline of the paper is as follows. In §2 we derive our main mathematical results for the 
amplification of an elliptical source. In §3, we derive analytical expressions for the entire lightcurve 
of a low mass (e.g. planetary) lensing event (§3.1) and for the wings of the lightcurve in general 
(§3.2). Sections 4 and 5 are dedicated to specific astrophysical sources. In §4 we calculate sample 
lightcurves for a giant star with a realistic limb-darkened profile and for an inclined accretion 
disk. In §5 we describe the time dependence of emission line profiles from a Keplerian disk during 
a microlensing event. Finally, §6 summarizes the main conclusions of this work. 



2. The Amplification Integral for Elliptical Sources 

In analysing the lensing geometry we use the two-dimensional sky coordinates. In this 
coordinate system, all projected distances are angular. 

We consider a general elliptical source with eccentricity e and semimajor axis a, which we set 
as our distance unit. The location of the lens at a given time can be defined by its displacement 
from the source center tq and by the angle «o between the source-lens vector and the major axis of 
the source. Figure 1 depicts the geometry associated with our notations. A lens with an Einstein 
radius e amplifies the flux from a point source at a distance r by a factor 

r 2 + 2e 2 

Mr) = r^—F ^ ■ (1) 

tVt 2 + 4e 2 

In the limit of geometric optics, an extended source can be viewed as a sum of infinitesimal 
point -like sources. Its amplified flux can therefore be obtained by weighting Aq(t) with the surface 
brightness distribution B(f) and integrating over the projected surface area of the source £5, 



F=a J m 7vwrw di: ■ (2) 

The net amplification is the ratio between this flux and the flux of the source in the absence of the 
lens, 

A = F/F , (3) 



where 



F = a 2 J B{r) d£ . (4) 



It is most convenient to describe the surface brightness distribution of the source in terms of a 
coordinate system which is centered on the source. We set the x-axis along the major axis of the 
ellipse (see Fig. 1), and adopt the following coordinates: x = pcosip,y = p\/\ — e 2 sin^. Curves of 
constant p are nested ellipses with semimajor axes p and equal eccentricity e ; p = 1 corresponds 
to the limb of the source. Note that if), which runs from to 2ir, is not the polar angle in the case 
of nonzero eccentricity. The area element is d£ = \A — e^pdpdifi, and so equation (Q) becomes 

2tt 1 

F = aVl-e 2 J f B(p,ip) palp dip . (5) 



In calculating the lensed flux in equation (||), it is advantageous to use lens-centered polar 
coordinates for which the amplification factor preserves its simple form ([l]). If we erect the 
corresponding x\ axis from the lens towards the source center, we obtain the coordinates: 
x% = Tcoscp, yx = Tsinc/) (see Fig. 1). To express the source boundary and brightness distribution 



in terms of t and (p, we need to transform the source-centered coordinates to the lens-centered 
coordinates. As seen from Figure 1, 

( x \ _ ( - cos a sin a \( X\ - t \ . . 

\ y I \ sina cosqq / I y\ I ' 



or equivalently, 



pcosip = — TCOs((j) + Oq) + Tq COS«o 



pyl — e 2 sin^ = rsin(0 + ao) — rosinao . (7) 

From these equations it is straightforward to express both p and ip in terms of r and (j). For p we 
get 

p = y 'u;i(0)t 2 - £j 2 (<£)t + a 2 , , (8) 



with 



p 2 

e ■ ,2/ 



wi(<P) •+- ^-sin + ao) 
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W2(^>) = 2t o [cos0+- ^sinQ sin((/) + a )l (9) 
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a = t W1+ - 2 sin^ao . 

The parameter uo is the value of p at the position of the lens (r = 0). If the lens lies within 
the source boundary then ctq < 1, and if it lies beyond the boundary, ctq > 1. For a given angle <f>, 
the integration limits in the r direction are found by setting p = 1 in equation (0), 



CJ 2 ± A/Wo + 40>l(l - On) 

r±(0 = 2uJl • ( 10 ) 

A line originating from the lens at an angle <f> enters the source boundary at r_ ((f)) (provided 
o"o > 1) and exits at t+ (</>). If ctq < 1, (j) runs from —ir to ir. To determine the (j) integration limits 
in the uq > 1 case, we set r + (0ii m ) = T_(</>i; m ) in equation (|T~0J) . After some algebra we find, 

/ COs20 Um \ , 4 , 4 2 2 ^ ^11\ 

Ui^AtaJ' 7 " - 2 ^»™ 2 °») (») 

_ (tS + <? - 2) ( T » " e . 2 <f 2a » ) ± 2^ + 7 - £j cos* a - 1 [ f^** )■ 

y e sm 2qo I v \ r o e cos 2"o / 

This expression generally yields four different solutions for 4>\\ m , due to the implicit ±7r ambiguity 
in this angle. The two values of interest are directed towards the source, and for them ui2{4>\im) > 0. 
We use 4>- G ( — 7r, 0) as the lower limit, and </>+ G (0, 7r) as the upper limit. Since the area element 



in the lens-centered coordinates is dH = rdr d(j>, the divergence of the integrand in equation 
at r = is eliminated. The explicit form of equation (12]) as a function of the lens position is: 



T+(4>) 



F(T ,a ) = a 2 J J BM)-j= 



r 2 + 2e 2 



-7T 



2 +4e 2 



dr 



and 






for ctq < 1 



for (Tq > 1. 



(12) 



(13) 



4>-T_(4>) 



So far we made no assumptions about the brightness distribution B(p,vfj). In the following 
we concentrate on elliptical source profiles with B being independent of ip. In particular, if the 
brightness can be written as a series in powers of p 2 



B(p) = Y l P i p 



2 1 



(14) 



i=0 



then equation (||) can be used to transform this series to a polynomial in r, 

k k—L Li \N 

b{ta) = Y,Y,^t?^- b{l+n) - 2L+n ■ 



L=0N=0 



L\ N\ 



(15) 



Here the constants B^ n > are the n-th derivatives of the brightness profile (|14|) with respect to p 2 
at the position of the lens, 
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P=o~0 i= n 



(i — n)! 



(16) 



Note that in this expression the derivatives of the polynomial (|l4|) are formally evaluated at the 
lens position even if it lies beyond the source boundary. 

By substituting expression (|l5|) into equations (|i"2|) and (fH|), we are left with r-integrals of 
the general form 



where /i = 0, ..., 2/c. These integrals can be solved analytically. We get 



Jq{p,v) = t\Jt 2 + 4e 2 



for p > 



h v +\(V;V) 



1 



2p + 3 



•2px Z^ 
.p/ i=0 



^^ + ^ 2 {2r 2 ^ 2 -^j:^)(-W e 2 r 



i T 2i 



(17) 



(18) 



(19) 



and for p > 



J2p(n,v)= ^pr 



t 2p+Vt2 + 4e 2 + Ap{ 2 p(-e 2 y +1 {arcsh£ + yV 2 + 4e 2 £ ^y 1 ^" 1 ! 



(20) 



In this way the two-dimensional integration over the source is reduced to a ^-integral. 
Equations (111) and (O) now read 



F(r ,ao) 4/EE ^^W]% (^) j2L+w( o, r+ (^)) # , for a < 1, (21) 



EE 

L=0N=0 

and 



^(ro,«o) = y / E E ^W);" 2 ^ iV i ? (^)j 2L+jV ( T _( ), r+W ) # , for a > 1. (22) 

Z J_L=0N=0 Ulyl 

For completeness, we substitute equation (|l4|) into equation (^) to get the unlensed flux 

F = ira 2 VT^J2^7 ■ ( 23 ) 

Surface brightness profiles of astrophysical sources are often defined through a set of discrete 
points [pi, Bi\, i = 1, . . . , n, with p\ = and /? n = 1. A straightforward application of the above 



results could be obtained through a least-squares fit of the form (14) to the data points. However, 
a satisfactory global fit for p £ (0,1) often requires the use of many high-order terms in the series. 
For the k + 1 terms of equation (|L1|) it is necessary to evaluate 2/c + 1 functions Jo • • • J2k> with 
increasing complexity at increasing k. As a result, an appropriate numerical integration scheme 
(such as Gaussian quadrature) would be more effective for high k. An alternative possibility is to 
interpolate between pairs of neighboring data points, so that for p G {pi, pi+\) the brightness is 
described by an adequate function Bi{p). A simple example for an interpolation scheme involves 

B i {p) = a l + b lP 2 , (24) 

where a^ and bi are determined by the continuity conditions at the boundary points, 

B l (p l ) = B l , Bi(p i+1 ) = B i+1 , i = l...n-l. (25) 

This interpolation scheme for the brightness profile does not have a continuous derivative at the 



data points pi. Smoothness can be easily achieved by adding a p 4 term to equation (24) and by 
adding conditions for the continuity of the derivative in equation (^5|) . 



Following the interpolation scheme of equation (|24]), the source is divided into n— 1 concentric 
bands, bounded by ellipses with semimajor axes pi. The brightness profile inside the j'-th band is 



described by values aj and bj. A line of constant angle <f) from the lens passes through a sequence 
of I bands in the r-direction, and crosses their boundaries at values n < r 2 < • • • < r^ +1 (see Fig. 



2). We correspondingly substitute in the integrand of equation (21) 

B( L + N U 2L+N (0,T + )-^J2 Bi s- +N)j ^+N(T u T i+1 ) , (26) 

where n = 0, t/ + i = t + , and Si is the index of the band crossed between Tj and Ti+i. An analogous 



substitution with ri = r_ should be performed in the integrand of equation (|22j). 



The sequence of bands that are crossed can easily be determined from the geometry illustrated 
in Figure 2. If the lens lies outside the source boundary (ctq > 1), one needs to find the innermost 
band which is crossed. The p-interval of this band contains the minimum value of p which is 
reached along the <p = const ray. Equation (||) yields 



Pmin = J O"o _ |77- • ( 27 ) 

If the lens lies inside the source boundary {<jq < 1), the sequence of bands starts at the lens 
position (p = do), and the band with p m in has to be determined only if uj 2 > 0. To calculate the 
crossing points Tj, it is then sufficient to know that intersections with ellipse pj occur at 



t= ¥ 1 . (28) 

L<jJ\ 

If the necessary number of bands is high, numerical evaluation of the r integral should be 
considered as an alternative. 

The method has so far been presented for an elliptical source. Most of the equations 
are considerably simpler for a circular source. By setting the eccentricity to zero, e = 0, the 
source-centered coordinates p, tp become regular polar coordinates, with a being the source radius. 
Equations @ reduce to: uo\ = 1, u>2(4>) = 2tocos0, and Co = To. The r integration limits in 
equation (|l0l) are 



t± = To cos(/> ± \j\ — Tq sin (j) ) (29) 

and the angular integration limits in equation ( |TT| ) are given by 

1 



± = ± arcsin — . (30) 

TO 



Due to the circular symmetry, the angular integration in equations ( |2l| ) and (22) can be performed 
over half the interval and the result multiplied by a factor of two. 

In principle, our method can also be used to describe the effect of spots on lightcurves of 
extended sources (Sasselov 1996). Since the observed flux is linearly dependent on the source 
brightness, one could simply superimpose the lightcurve of an elliptical spot on top of the 
background lightcurve. 



9 



3. Analytical Results 

3.1. Low Mass Lenses 

The results obtained for brightness profile (|14) can be used to study the limiting case when 
the Einstein radius of the lens is much smaller than the source size (e <C 1). The lensed flux in 



equations (|2l| ) and (22) depends on e only through the functions Jh in equations (|18|)-([2C|). For 
small e, these functions can be expanded as power series in e. We make use of the expansions 



Vt 2 + 4e 2 



r + 



2e 2 



2e 4 



+ 0(e 6 



arcsh- 



2e 



lne + lnr + ^ + 0(e 4 



If the lens position is inside the source boundary (o"o < 1), we obtain 



L > 2 



J (0,r + 
Ji(0,r+ 

•7 2 (0,7+ 



2e 2 



2f 4 



2 o 8e 3 

3 r + + lT 






+ ^ 4 ) 



T\ 



4e 4 lne + (41nr+-3)e 4 + o(e 4 



_L+2 



4e 4 



L + 2 



L 



r^ + (6 4 ) 



If the lens is outside the source (o"o > 1), we get 

2 



L/2 Jl(t_,t + ) 
J 2 (t-,t + ) 



T i+2 , 4e L-2 

L+2 L-2 

4 l r + 



T + 



+ 0(6 4 ) 



r 



+ 4e 4 lnr 



+ o(6 4 ) 



The observed flux can now be obtained from equations (21) and (22): 

F(r ,a ) = ^o + 2vra 2 e 2 L ,(0) + 0(e 4 lne) for a < 1, 

and 



F(t , oo) = Lo + 2e 4 a 2 J^dE + 0(e> 



for o"o > 1. 



(31) 



(32) 



(33) 

(34) 
(35) 



Since L?' ) is the brightness at the position of the lens, the lightcurve traces the brightness 
profile along the lens trajectory to leading order (e 2 ). For gq < 1, the excess flux equals twice 
the Einstein ring area times B(°\ The accuracy of this leading-order approximation will be 
demonstrated through a numerical example in §4. We expect the first two terms in equations (531 
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and ( p5D to hold in general for an arbitrary surface brightness distribution, as long as the variations 
in B(r) are weak on the scale of the Einstein radius of the lens. Note that the convergence of the 
above expansions requires that 2e be smaller than t_ and r + . Hence, the derived expression for 
the lightcurve is valid as long as the lens lies a distance greater than 2e away from the source 
boundary. 



3.2. Lightcurve Wings 

We next study the effect that a finite source size has on the wings of a general lightcurve. For 
this purpose we start from equation (0) and place no restrictions on the brightness distribution 
B(r). We associate the "wings" of the lightcurve with the region over which the lens-source 
separation is larger than both the Einstein diameter of the lens and the size of the source, i.e. we 
assume that t 3> 2e and t 3> 1. Rewriting equation ([|) in terms of source coordinates, we get 

2e 2 \ / 4e 2X 2 








F = aVl - e 2 / / B(p, </,)! + _ 1 + — p dp dip , (36) 



where r 2 can be expressed from equations (^), 



t = r — 2 / 0To(cosao cos ip — vl — e 2 sin ao sin ip) + p (1 — e sin ip) . (37) 

The amplification factor in the integrand of equation (^) can be expanded in powers of Tq , 




2 2g 4 8oe 4 

= 1 H j- H r- (cos «o cos -0 — vl — e 2 sin ao sin ^) + 0(tq 6 ) . (38) 

r o r o 

By substituting this result into equation ( |3"6| ) we obtain, 

/ 2e 4 \ 8e 4 a 2 Vl — e 2 /" /" 

F = I 1 H j- Fq H g B (cos ao cos ^ — vl — e 2 sinaosin'0)p d^ dp + ©(tq - ). 

V / o Q 

(39) 

The leading correction to the unlensed flux (of order Tq ) is independent of the source 
structure, and appears also for a point source. The next term (of order t ~ ) vanishes for a 
point source, as well as for any symmetric source with B(p,tp) = B(p). The integral in this 
term is, in fact, the projection of the "brightness dipole moment" / rB{r) dT, on the source-lens 
direction. This term therefore reflects the asymmetry in the brightness distribution of the source, 
which could result from the existence of hot or cold spots on its surface (Sasselov 1996; Gould 
& Miralda-Escude 1996). The term vanishes if we redefine the positions relative to the center of 
brightness, rather than the geometrical center of the source. 
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4. Sample Lightcurves of Astrophysical Sources 

The microlensing lightcurve of the source is obtained from the time dependence of the lens 
position coordinates, tq and ao- A linear motion of the lens with respect to the source can be 
described by two parameters: the angle /3 between the major axis of the source and the lens 
trajectory, and the impact parameter p (see Fig. 3). To avoid ambiguity, we assign a negative sign 
to p if the source center lies to the right of the lens trajectory. For convenience, we normalize the 
time t in units of semimajor axis crossing-time and set t = at closest approach. The geometry 
illustrated in Figure 3 then yields 



to = ^P 2 + t 2 (40) 

cosao i 1 / cos/3 sin/3 



sin ao I tq \ — sin /3 cos /3 




(41) 



In the case of a circular source (e.g. a lensed star), only p needs to be specified and the angles (5 
and ao are redundant. 

As an example of stellar microlensing we present a theoretical lightcurve for the MACHO 
95-30 event (Pratt 1996, Alcock et al. 1997). The source, a red giant in the galactic bulge, was 
amplified by a lens with an Einstein radius e ~ 13 and an impact parameter p ~ 0.7. The R-band 
brightness profile shown in Figure 4 is based on an atmosphere model calculation of the source 
by Sasselov (1996). The computed lightcurve is plotted in Figure 5 together with a point-source 
lightcurve, for comparison. 

Figure 6 compares the exact lightcurve of a low mass lens with the leading order term in the 
expansion described by equation (p4|). In difference from the real MACHO 95-30 event, we assume 
in this example that the lens has an Einstein radius e = 0.1 and an impact parameter p = 0.2. In 



this case, the leading order term (of order e 2 ) in equation (34) captures the main characteristics of 
the lightcurve to a reasonable accuracy. 

The case of an elliptical source can be illustrated through the example of microlensing of an 
AGN disk. A circular disk which is inclined by an angle i relative to the line-of-sight, would appear 
as an ellipse with eccentricity e = sini on the sky projection. The semimajor axis would then 
equal the disk radius, and the p = const contours would be elliptical projections of circles, with i/j 
being the standard polar angle on the face of the disk. In calculating a sample lightcurve for this 
application, we use the brightness profile of a stationary thin accretion disk (see Frank, King, & 



Raine 1985) with B(p) oc 1 — \/pi n /p\ p ■ We assume a disk with an inner radius pi n = 0.1 and 
an inclination i = 60°, which is lensed by a lens with an Einstein radius e = 5. Figure 7 contains 
lightcurves for events with the same impact parameter p = 0.8 but different approach angles (3. 
While the lightcurves of a circularly symmetric source are invariant to /?, the inclination of the 
disk leaves a distinct /3-dependent signature on the amplitude and symmetry of its lightcurve. 
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5. Microlensing of an Emission Line From a Keplerian Disk 

Next we consider the effect of microlensing on the spectral shape of an emission line from an 
accretion disk. Because of the Doppler shift associated with the disk rotation, the amplification as 
a function of wavelength could help unravel the inclination and surface brightness distribution of 
the disk. 

To examine the microlensing distortion of the line profile we adopt the simplest model of a 
planar thin Keplerian disk with an inner radius pi n and inclination i. We parametrize the disk 
surface by the coordinates p, ip which were introduced in §2, and assume a Keplerian velocity 
profile v(p) = Vi n \/pin/p, where Vi n is the velocity at the inner edge of the disk. We assume that 
the cold gas in the disk emits a line at a rest-frame wavelength Ao, and with a thermal width 
that is much smaller than the Doppler width induced by the rotation of the disk. We therefore 
approximate the rest-frame line profile by a delta-function centered on Ao- For simplicity, we 
neglect any relativistic effects on the radiation emitted by the disk. Due to variations of the 
line-of-sight velocity across the source, the Doppler-shifted wavelength reaching the observer from 
a point (p, ip) is 



A = A r 



v(p) . . ,' 
1 H sin i cos ip 



(42) 



where c is the speed of light. In this notation, ip = is on the redshifted half of the disk. The 
spectral intensity profile of the line emitted by a microlensed disk is then 



1(A) = j(p)5 



\ \ / i j U W • • / 
A - Ao H sin i cos ip 



Aq(t(p, ip)) pcosidipdp , (43) 



Pin 



where j(p) is the spectral emissivity and Aq(t) is the amplification factor in equation ([!]). The 
displacement r(p, ip) of the lens from the integration point is given by equation (37) with e = sini. 
We normalize the wavelength shift by its maximum value, (vi n /c) sini, and denote the normalized 
fractional shift by A = (A/Ao — 1) /[(^m/ c ) sini]. Integration of the delta-function in equation ( |43"D 
over ip yields 



pmax 
Vr 



A (T(p,aiccoB[Ay/p/pi n ]j) + A ( r(p, - arccos[A^p/ p in })) 

I A [l + A^sini] = / j(p)—i — ^=7=r3f I ■ L Pdp- 

V c J J M{Vin/c) VPin/p- A tarn 

Pin 

(44) 
The two terms in the numerator reflect the fact that points with angles ip and —ip have the same 
line-of-sight velocity. The upper integration limit depends on wavelength. For small wavelength 
shifts |A| < y/pi n , one gets contributions from all radii so that p m ax = 1- Larger shifts |A| > y/pi n 
originate only from the inner disk, with p ma x = Pin I A 2 . 

For concreteness, we assume a power-law profile for the disk emissivity, j(p) oc p n . Figure 8 
shows the sensitivity of the unlensed line profile to the power-law index n. For the purpose of this 
illustration, we have kept the total luminosity of the disks constant, namely J j{p)pdp = const. 
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The unlensed profiles show two characteristic peaks at |A| = yjpi n . We have chosen pi n = 0.01 
and i = 60° in the examples shown. As is apparent, the more centrally concentrated the emissivity 
law is, the more pronounced the line wings are and the less pronounced the peaks are. This 
tendency results from the enhanced weight that the inner region of the disk obtains in the intensity 
integral (44), as n gets more negative. The wings of the line correspond to high rotation speeds 



which are generically obtained in the inner region of the disk. 

Figure 9 shows the effect of lensing on the line profile for n = — 1. In this calculation we 
assumed a lens with an Einstein radius e = 1, located at several different positions along the major 
axis of the source (ao = 0). For clarity, we show results only for the redshifted side of the disk. 
Apart from an overall amplification of the total flux in the line, lensing introduces an asymmetry 
between the intensities of the two peaks. The position of the peaks remains unaffected. The curve 
with ro = 0.1 illustrates another feature, present when the lens is positioned on top of the disk. In 
this situation, the profile diverges at the emission wavelength of the geometric point behind the 
lens, A = \J Pm I Picns cos V'lens- This divergence is integrable, and so a filter with a finite wavelength 
resolution would register it as a third peak. In reality, the formal divergence will be smoothed 
out by the thermal width of the line. The time-dependent intensity profile around the moving 
third peak can be used to map the axisymmetric brightness profile of the disk as the lens scans 
the source along its linear trajectory. 

The asymmetry of the microlensed line profiles in Figure 9 indicates a change in the mean 

redshift of the emitted line, depending on the lens position. This effect has been previously 

studied for various models of the broad line regions of quasars (Nemiroff 1988; Schneider & 

Wambsganss 1990). The redshift change Az can be simply computed from the spectral intensity / 

in equation (|44]) as 

l 

/ IAdA 

Az = — sini . (45) 

JIdA ° 

-i 

Figure 10 shows the change in redshift as a function of lens position along the major axis of the 
disk, for the same situation as in Figure 9. In this case the redshift can change by more then 10% 
of the maximum Doppler shift from the inner edge of the disk, if the lens is positioned close to the 
inner edge. 



6. Discussion 

We have developed an efficient method for computing microlensing lightcurves of sources 
with an elliptical symmetry and a general surface brightness distribution. Equations (|l8|)-(f23"|) 
express the amplification of elliptical sources in terms of a one-dimensional integral. This integral 
is considerably simpler for a circular source. The method is well suited for the study of spectral 
changes due to microlensing. Such changes have already been observed in the MACHO 95-30 
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event (Alcock et al. 1997). We are currently using the technique for modelling these observations 
(Sasselov, Heyrovsky & Loeb 1997). The method can also be easily extended to describe the 
effect of elliptical spots on the background lightcurve.We have derived fully analytical results in 



the limit of low lens masses [cf. equations (p3), (35)] and for the wings of the lightcurve in general 



[equation (p9|)]. In both cases, the lightcurve carries information about deviations of the source 
from elliptical symmetry (e.g., due to hot or cold spots). The low mass limit case applies to lensing 
of giant stars in the Milky Way bulge by planets of mass M p \, for which the Einstein radius is 
~ 3i?0 x (Mpi/M^) 1 ' 2 . The lightcurve in this case obtains the form of a low amplitude increase 
in the source intensity over the source crossing time. The amplitude of the excess flux is twice the 
area of the Einstein ring times the brightness of the source at the lens position [cf. equation ([34}) ]. 
Events of this type might soon be discovered through ongoing microlensing searches for planets 
(Albrow et al. 1996; Pratt et al. 1996). Finite source signatures appear whenever the impact 
parameter is comparable to the source size, including high amplification events caused by massive 
lenses. Figure 5 shows the theoretical lightcurve corresponding to the source and lens parameters 
of a microlensing event of this type, MACHO 95-30 (Pratt 1996, Alcock et al. 1997). The 
detection rate of events with finite source signatures is expected to increase if future searches will 
focus on the source population of giant stars (Gould 1995a). When such signatures are combined 
with spectroscopic identification of the sources, they can be used to break the degeneracy of 
microlensing lightcurves (Loeb k, Sasselov 1995; Gould & Welch 1996). 

Our formalism is particularly suitable for the analysis of microlensing of inclined accretion 
disks which appear elliptical in projection on the sky. The microlensing lightcurve of the disk 
depends on its surface brightness distribution and inclination (cf. Fig. 7). If a thin disk emits a 
spectral line, then it is possible to deconvolve its properties from the variation of the line profile 
during a microlensing event. When the lens position overlaps with the disk, a third peak is added 
to the standard "double-peak" profile of the disk emission. The microlensing peak appears at the 
emission wavelength corresponding to the Doppler velocity behind the lens (cf. Fig. 9). As the 
lens moves along a straight line, the temporal evolution of the intensity and wavelength of the 
third peak can be used to map the disk structure. Thus, microlensing of a spectral line offers a 
unique opportunity for unraveling the surface brightness and velocity profiles of accretion disks. 
The main practical obstacle for such mappings is the potential for intrinsic variability in the 
disk properties during the event. Figure 10 illustrates the change in the mean redshift of the 
spectral line, due to the asymmetry of the microlensed profiles. The change can exceed 10% of the 
maximum Doppler shift in the disk. 

The signatures described in this paper should be more complicated in environments where the 
microlensing optical depth is not small, such as the cores of galactic macro-lenses at cosmological 
distances. However, our isolated point lens calculation still applies to the outer halos of such 
galaxies, where the optical depth is low. Quasar absorption spectra can be used to select such 
foreground galactic halos. In particular, damped Lya or metal-line absorption features are often 
associated with projected impact parameters of ~ 10-100 kpc from galactic centers (Steidel et al. 
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1994, 1995, 1996), and could therefore be used to signal an intervening galactic halo. If galactic 
halos are composed of Massive Compact Halo Objects, then a significant fraction of all damped 
Lya absorbers should show evidence for isolated microlensing events (Perna & Loeb 1997), of the 
type discussed in this work. Isolated microlensing events from intergalactic stars are also possible 
but only for a minority of all quasars (Dalcanton et al. 1994; Gould 1995b). 

We thank Dimitar Sasselov for providing the brightness profile of the MACHO 95-30 
source, and Robert Nemiroff for suggesting the calculation of the change in the mean redshift of 
microlensed emission lines. We also thank both for valuable comments on the manuscript. This 
work was supported in part by the NASA ATP grant NAG5-3085 and the Harvard Milton fund 
(for AL). 
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Lens 



Fig. 1. — Geometry of the lens and an elliptical source projected on the sky. r marks an arbitrary 
point on the source. 
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Fig. 2. — Sketch of the integration method geometry. The source is divided into three bands, each 
of which is characterized by its separate section of the brightness profile interpolation. For the 
angle 4> m this figure, the band sequence is {s±, S2, S3} = {3, 2, 3} and t\ = , T4 = r+(0). 
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Lens 



Fig. 3. — Trajectory of a lens in linear motion with respect to the source. 
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Fig. 4. — Brightness profile of the MACHO 95-30 source in the R-band, normalized by its central 
value. The profile is based on a model atmosphere calculated by Sasselov (1996). 
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Fig. 5. — Simulated lightcurve of the MACHO 95-30 event (solid line), using the brightness profile 
in Figure 4. The lens was chosen to have an Einstein radius e ~ 13 and an impact parameter 
p ~ 0.7. For comparison, the dashed curve shows the point-source lightcurve for the same event 
parameters. 
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Fig. 6. — Lightcurve of the MACHO 95-30 source for a lens with an Einstein radius e = 0.1 and 
an impact parameter p = 0.2. The exact lightcurve (solid curve) is compared to its leading order 
approximation in the low lensing-mass limit (dashed curve). The brightness profile of the source 
was taken from Figure 4. 
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Fig. 7. — Lightcurves for an accretion disk inclined by an angle i = 60°, which is microlensed by a 
lens with an Einstein radius e = 5. The three curves correspond to events with the same impact 
parameter p = 0.8, but with a different lens orientation. The solid line corresponds to a trajectory 
parallel to the minor axis of the source {(3 = 90°), the dotted line to a trajectory diagonal to the 
axes {(3 = 135°), and the dashed line to a trajectory parallel to the major axis {(3 = 180°). Zero 
time corresponds to closest approach in all cases. 
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Fig. 8. — Dependence of the unlensed profile of a spectral line emitted by a Keplerian disk, on the 
radial power-law index of the disk emissivity. The different line profiles correspond to power-law 
indices of n = — 1 (solid line), n = —2 (dashed line) and n = —3 (dotted line), for a disk with an 
inner radius pi n = 0.01 and inclination i = 60°. The vertical axis is in arbitrary units. 
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Fig. 9. — Effect of microlensing on the line profile for a disk with n = — 1 in Figure 8, and for 
various lens positions along the major axis of the projected disk. The lens has an Einstein radius 
e = 1, and is located at to = (dashed line), to = 0.1 (dot-dashed), and to = 1.2 (dotted). The 
unlensed profile is plotted for comparison (solid line). 
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Fig. 10. — Change in mean redshift of the emitted line as a function of lens position tq along the 
major axis of the projected disk. The change in redshift is plotted in units of maximum Doppler 
shift from the inner edge of the unlensed disk. Disk parameters are the same as in Figures 8 and 
9. The three marked points correspond to the lens positions for the line profiles in Figure 9. 



